On the Harmonic approximation for large Josephson junction coupling charge qubits 
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We revisit the harmonic approximation (HA) for a large Josephson junction interacting with some 
charge qubits through the variational approach for the quantum dynamics of the junction-qubit 
coupling system. By making use of numerical calculation and analytical treatment, the conditions 
under which HA works well can be precisely presented to control the parameters implementing the 
two-qubit quantum logical gate through the couplings to the large junction with harmonic oscillator 
(HO) Hamiltonian. 
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£L|' I. INTRODUCTION 

<D 

For quantum information science and technologies, it is crucial to build the fundamental quantum logic gates 
Together with the basic single bit logic gates, the non-trivial two bit gates constitutes the fundamental blocks for the 
quantum network of quantum computing . Physically, the two-bit quantum logic gate are based on the entanglement 
of two-qubit through the controllable interaction between them [3J- 

In most scenarios to build the two-qubit logic gates, two qubits interact with a common object as data bus. When 
the normal frequency of the data bus is off resonance with respect to the energy level spacing of the qubits, the 
variables of the data bus can be removed and then induce an effective interaction between the two qubits, which can 
create an entanglement based on two-qubit logic gate. Most recently, this idea was used to construct the two-qubit 
logic gates of the superconductor Josephson junction where the data bus is implemented as a large junction. Since 
the recent experiments Q have shown longer time quantum coherence, the Josephson junction qubits ( JJ qubits ) 
including charge qubits, flux qubit and the single junction qubit, the practical schemes to implement the two-gate 
operations of J J qubit become very important. We notice that many of these schemes depend on the approximate 
approach to treat the large junction as a harmonic oscillation. We call this approximation the harmonic approximation 
(HA). Due to the linearity of its coupling to the two charge qubits, the variable of the harmonic oscillator (HO) is 
easy to remove in a dynamic way or adiabatically [(j. 

Now the validity of HA becomes the focus in search of the implementation of the JJ qubit based quantum compu- 
tation. In this paper, we will tackle this problem by the numerical calculation and the analytical considerations. 

cr 



o 

^— > ■ 



II. THE MODEL FOR THE COUPLING OF A LARGE JUNCTION TO A CHARGE QUBIT 



We consider a system of coupled charge-phase qubits 0, H| as illustrated by the electronic circuit shown in Fig. 1 . 
We are interested in the limit case that the large junction (of capacitance C) stays at a low excitation or even a 
thermal state. The Cooper pair box is a small junction of capacitance C (C = C") that form a superconducting 
loop. C g is the capacitance of the gate, and Ej is the Josephson coupling energy. The total Hamiltonian containing 
the reduced Coulomb energy and the three Josephson coupling energies reads 

H = E c (n -n g ) 2 - Ej (cos (p + cos ip") + E c n 2 - Ej cos 9 (1) 

where E c — 2e 2 /Cs, n g = C g V g /2e, E c — 2e 2 /C and Cs = C g + C + C". Here n is the number of Cooper pair 
on the island, while n is the number of Cooper pair on the Coulomb island connected with the large junction, ip , 
ip" , and 9 are superconducting-phase differences across the relevant junction. They are related through the fiuxoid 
quantization condition around the loop 9 + ip — ip" — 20 and 0/. = Ti^> x /(j)Q. Introducting <p = (ip + ip")/2, we rewrite 
the Hamiltonian as 

' ' 

H = E c (n - n g ) 2 + E c n 2 - Ej cos 9 - 2Ej cos(- - 9) cos ip (2) 
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FIG. 1: Schematic electric circuit diagram of two charge-phase qubits coupled through a large junction. 
For charge-phase qubit, when n g — 0.5 and Ej is not much larger than E c , the avove Hamiltonian can be reduced 



H = E c (n 2 - T cos 8 + T sin - • cr 3 ) 



(3) 



where T — Ej/E c , T = E'j/E c and is Pauli matrix. The operators n and 9 satisfy the quantization condition 
[9,n]=i. 

To consider the validity of the HA we only focus on the large junction together with its coupling to the charge 
qubit. In this sense we write down the corresponding Hamiltonian 



H = n 2 - T cos 6 + T 1 sin -a 3 

#(+) 




(4) 



bmce the relation H^+\9) = H^(-6) holds, we know that these two Hamiltonians H^ + \H^ > have the same 
spectrum. Thus one can get all the eigenstates of (1) from the eigenstates of H^ + \ In the following, we will 
concentrate on the eigenstates of the Hamiltonian in the qubits states spanned by discrete charge states |n >, 
n = 0, ±g, ±1, ±|, ±2, • • • , ±oo. The periodic potential can be expanded around the equilibrium points 6o, where 
sin(#o/2) = ±T'/4T at which potential has minima. Then can be written as 



^2 , 16T -T ' 



+2arcsin(-T74T))] 2 -^A 
6 G [(2i - 1)tt, (2i + l)ir],i = even 



(•5) 
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32T 
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-2arcsin(-T'/4T))] ; 



&T A +T 
8T 



I 9 e [(2i - 1)tt, (2i + 1)tt], i = odd. 



In large T limit, the tunnelling between two neighbor wells is forbidden. The eigen functions of ajai are standard 
Harmonic wave functions with the argument [9 — (2iri ± 2 arcsin(— T 1 /AT))] when i is even or odd. Obviously, the 
spectrum of the Hamiltonian [jj is the same as the simple harmonic oscillator, i.e. 
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16T 2 - T' 2 T' 2 
8T ~ 8T 
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0,1,2, ), 



3 



but the eigcnstates are multi-fold degenerate. Notice that the above discussion is based on the assumption that 
there is no tunnelling between two individual parabolic potential wells. The validity of this assumption need to be 
investigated analytically and numerically. 



III. VARIATIONAL SOLUTION FOR THE MOTION EQUATION 

In this section, we solve the Schrodinger equation governed by large junction Hamiltonian (2). In the Hilbert space 
spanned by the eigen state \p > of n, where n\p >= p\p > and p£ [—00,00] the Hamiltonian can be written as 

/oc 
\p X p\p 2 dp 
-OO 

rp POO 

-77 {\p+l><p\ + \p><p+l\)dp (7) 

^ J-00 

T' f°° 1 1 
-*~2 J (\P+2 ><P\-\P><P+7 l \)dp. 

Based on the transformation p — > n + k, dp — ► ^ n j^ 2 dk, where n is integer and half-integer, this Hamiltonian 
is written as 



H = i / II '' (//, 



1/2 



1 f ' 

= ~2J / [\n + k >< n + k\(n + k) 2 
„ Jo 



T 

--(\n + k+l><n + k\ (8) 

+ \n + k ><n + fc + l|) 

T 1 
-i—{\n + k+ - ><n + k\ 

-\n + k >< n + k + ^\)]dk. 

The Hamiltonian is invariant under the transformation k — > N/2 + k, where N is integer. It indicates that the whole 
space can be decomposed into the invariant subspaces denoted by k. Therefore the eigenstates of (4) can be written 

as 

>= n \^ k> ( 9 ) 

0<fe<l/2 

This means that one can obtain the eigenstates in any subspace from \tp > by the projection operator Pk = J2 n \n + 
k X n + k\, whenever Pk\tp >^ 0. What we concern is the eigenstates in k = subspace. Since it is not so easy to 
tackle the eigen problem analytically in a subspace spanned by discrete basis, our strategy is that we first study the 
eigenstate of (3) by variational approach, then project the eigenstate onto the k = subspace. 

In order to investigate the eigenstates of the Hamiltonian in large T" limit, we introduce a set of trial harmonic 
wave function, 



[°° 1 

\<p n >= N n / exp(--a 2 p 2 + if3p)H n (ap)dp \p) 
J-00 z 



(10) 

where a, (3 are undetermined constants, N n = ^^ n \ 1S rcnormauza tion factor. We assume a -C 1. A direct 
calculation yields the following matrix elements 



(ip m \H\<p n > = AS m<n + B{ sj m(m + V)6 m , n +2 + \J n(n + l)(5 mi „_ 2 ) (11) 
+iC(m«— <V„+i - na——d m . n - 1 ) 
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where 

2m +1 ,ma 2 T I 



+ (^^ - l)T'sin^-exp(--f-a 2 ) (12) 



2a 2 ' v 2 ""^ ""^ A 

,ma 2 . /3 .1 

~2 exp( Te 

I a 2 „ „ ,1 



s = I^ 4Tcos/3exp( ~I a2) ( 13 ) 

+T'sin^exp(-^ a 2 )] (14) 

1 T" (3 1 

C = Tsin/3exp( — a 2 ) H cos — exp( a 2 ) 

4 2 2 y 16 

The trial wave functions are approximate wave functions of the Hamiltonian if all the off-diagonal elements vanish. 
Taking B = C = 0, we have 

. /3 T' 4 32T 

Sm 2 =4T' a = 16T 2 -T' 2 - (15) 

The spectrum is 



1 / 16T 2 - T' 2 T' 2 
Er n = (m + -) ] /— T--. (16) 

What we concern is the eigenstate in the subspace k = 0, which can be obtained by the project operator P = 

Enl"><»l 

^2 \n >< n\f m > 

n 

= iV m ^exp(-in 2 a 2 + i(3n)H m (an) \n) (17) 

n 

where (3 — 2 arcsin( ^ ) and 2n — 2 arcsin( ^ ) which correspond to the minima of the potential T cos 9 + T' sin | . Then 
in \n > representation, the lower approximate eigen wave functions are 



(18) 



±1 ±1,±|,±2,--- ,±oo. 

Notice that each energy level in (16) corresponds to a narrow energy band for the exact spectrum of the Hamiltonian 
(2). Employing the Hellmann-Feynman theorem, one can evaluate the band width 

/•1/2 CiTT 1 

AM = jf (^)mdk = - (19) 

which is much less than the gap between two neighbor bands. The eigenstates of H k share only one level from each 
band and the approximate wave function \<p n > is just the linear combination of all the eigenstates corresponding to 





= iV m exp(-^a' 
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= 0,1,2,- •• ,n = 



IV. NUMERICAL RESULTS 



In order to verify the validity of HA, numerical calculation is employed to compute the energy levels and wave 
functions. Since the whole Hilbcrt space is constitutive of invariant subspaces and in the practice system the charge 




FIG. 2: The probability amplitudes of three lowest levels for T — 100, (a) T' = and (b) T' = 60 obtained by numerical 
and variational methods. Solid ( open ) square, circle and triangle denote numerical ( variational ) results of eigenstates of 
n = 0,1,2. 
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FIG. 3: The energy levels for T = 100, (a) T' = and (b) T' = 60 obtained by numerical ( solid square ) and variational ( 
open square ) methods. It shows that the value of T' doesn't effect the results obtained by two methods so much which agree 
quite well for lower levels. 
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number n is restricted to be integer and half-integer, exact diagonalization is performed to get the lower energy 
spectrum and corresponding eigen functions of the Hamiltonian in k = subspace. In Fig. 2 and 3, the energy 
levels for systems with T — 100, T' = and T = 100, T" = 60 obtained by numerical and variational methods are 
plotted. It shows that the HA is valid for wide range of T 1 . The sudden deviations occur around zero energy since 
the tunnelling effect between neighbor wells is non-neglectablc. Under the threshold point numerical and variational 
results agree well. In Fig. 4 and 5, the probability amplitudes of three lowest levels for systems with T = 100, T" = 
and T — 100, T' = 60 obtained by numerical and variational methods are plotted. It shows that numerical and 
variational results agree well. 
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